Naturalness of CP Violation in the Standard Model 
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We construct a natural measure on the space of CKM matrices in the standard model, assuming 
the fermion mass matrices are randomly selected from a distribution which incorporates the observed 
quark mass hierarchy. This measure allows us to assess the likelihood of Jarlskog's CP violation 
parameter J taking its observed value J ~ 3 x 10 -5 . We find that the observed value, while well 
below the mathematically allowed maximum, is in fact typical once the observed quark masses are 
assumed. 



I. INTRODUCTION 

Kobayashi and Maskawa's beautiful explanation for 
why CP violation is naturally expected in the standard 
model of particle physics, with three or more families [l|, 
has been recently rewarded with a Nobel prize. Still, 
the magnitude of the observed CP violation remains un- 
explained, and this is a major task awaiting a funda- 
mental theory. Even before we have such a theory, it 
is interesting to ask whether the observed CP violation 
is "finely tuned," and therefore poses a major new puz- 
zle, or whether it is in some sense "typical," of what a 
fundamental framework might be expected to predict. 

According to Kobayashi and Maskawa, the observed 
CP violation originates in the Yukawa couplings between 
right- and left-handed quarks. These couplings are not 
entirely physical because they can be altered by a global 
redefinition of the basis of fermion fields. To resolve this 
ambiguity, Jarlskog has constructed a basis-independent 
measure of CP violation, called J. In this Letter, we 
construct a natural measure on the space of mass ma- 
trices, assuming the observed quark mass hierarchy but 
predicting the CKM mixing angles and phase. From this 
we infer the typically expected value for J. We find that 
the observed value of J, while far smaller than the max- 
imal mathematically allowed value, is indeed very likely. 
Hence we conclude that, at least with the measure we 
define, there is no "fine tuning" problem in this aspect of 
the standard model. 

We begin from the low energy effective Lagrangian de- 
scribing standard model quarks. In the gauge basis, this 
contains mass terms 



= -Q L mQ R - Q' L m' Q' R + h.c, (1) 

where Q = (t, c, u), Q' — (6, s, d) and m and ml are arbi- 
trary 3x3 complex matrices, with dimensions of mass. 
Following 0, we define the dimensionless mass matrices 
AI = m/A and M' = m'/A', where A and A' are mass 
scales which may be chosen for convenience. Via a re- 
definition of Qr and Q' Rl M and M' may be rendered 
Hermitian, with non-negative eigenvalues. They can be 
represented as follows 

M = U^DU, M' = U n D'U', (2) 



where D and D' are real and diagonal, D = 
diag(M t , M C ,M U ) and D' = diag(M 6 , M s , M d ). The fi- 
nal step is to redefine UQl = Q™ ass and similarly for 
Qr-iQl-iQri removing U and U from the mass terms. 



They now appear only in the charged current interaction 
terms in the Lagrangian, J£ c — Q L j^VQ'™ 3 * 5 



h.c, 



with V — UU'^ being the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix. 

The normalized mass eigenstates are only defined up 
to a phase. Jarlskog defined a measure of CP violation 
which is invariant under the choice of these phases [2], |3J] , 



J = -i det[M, M'] /(2TB) , 



where 



T = (M t - M u )(M t - M C )(M C 
B = (M h - M d )(M b - M S ){M S 

In terms of the CKM matrix V, 

J = Im (V n V 22 V{ 2 V 2 \) 



M„), 
-M d ). 



(3) 



(4) 



(5) 



and | J | has an elegant geometrical interpretation as twice 
the area of the so-called unitarity triangle. 

Mathematically, © represents the space of mass 
matrices M (with three distinct eigenvalues) as 
(£/(l) 2 x e 3 )\SU(3) x K 3 , for generic D, because left- 
multiplication of U by a diagonal phase matrix P 6 
SU(3), i.e., an element of U(l) 2 , does not alter M. 
To avoid overcounting, we must in addition factor out 
by elements of 5[/ (3) which permute the masses, as we 
shall detail below (S3 is the symmetric group on 3 el- 
ements). Because only the product UU'* appears in 
physical quantities, the two SU(3)'s are redundant and 
we may identify them with a diagonal SU(3) subgroup. 
Thus the final physical space of angles and phases en- 
coded in V may be thought of as the double coset 
(U(l) 2 x e 3 )\SU(3)/(U(l) 2 x 63). Our goal is to con- 
struct a natural measure on this space which incorporates 
the observed values of the quark masses. 

The standard measure on the space of Hermitian ma- 
trices, DM, will be detailed shortly. Since the space 
of eigenvalues is not compact, we must also include a 



2 



weight function which makes the integral over the eigen- 
values converge. We choose the weight function to distin- 
guish the three up- type quarks and give them appropriate 
expected masses, and similarly for the three down-type 
quarks. The simplest assumption is that the measure 
factorizes, 



DMf (Tr{M 2 A)) DM'f (Ty{(M') 2 A')) 



(6) 



where A and A' are Hermitian and positive definite, and 
we shall further assume [A, A'] = 0. 

Since physical quantities are invariant under U — > 
UW, U' — > U'W, without loss of generality we can choose 
A and A' to be diagonal. For simplicity and ease of tech- 
nical calculations, we shall choose the function / in © 
to be a decaying exponential so M and M' are governed 
by Gaussian distributions. Our proposal is to fit the di- 
agonal matrices A and A' to the observed quark masses, 
and then employ the resulting distribution given by ^ 
to infer a probability measure on J. 

The introduction of the matrices A and A' means that 
the invariance of the measure DM DM' under separate 
conjugation of M and M' by arbitrary elements of J7(3), 
i.e., under the action of 17(3) x U(3), is broken down to 
the action of the diagonal subgroup U(l) 2 x U(l) 2 which 
commutes with A and A'. We find that this symmetry 
breaking is necessary to obtain a distribution that re- 
produces different expectation values for squared quark 
masses. 

Full details of the calculations presented in this Letter, 
as well as a more general discussion of possible measures 
on the space of CKM matrices, may be found in the sep- 
arate publication 



II. CALCULATING THE MEASURE 

The natural line element on the space of Hermitian 
matrices is 



ds 2 = Tr(dMdM) 



(7) 



Tr 



(dDdD) + 2Tr((cff/C/ t D) 2 - (dUU^f 



D- 



which is invariant under conjugation under U(3). Our 
parametrization for the matrix U, an element of the space 
U(1) 2 \SU{3), is 



U = WT, 



R ■ 



(8) 



where 



W = 



and 



CyC Z 



CyS Z 



jx°y*-z ^x^z 



~x°y^z 



^x^y 



Sx^z 6 C x SyC z S X C Z 6 C x SyS z C x Cy 



T R = 



diag( 



gir+it gii — it g — 2ir 



(9) 



(10) 



where s x = sinx, c x = cosx, etc., and the ranges of the 
coordinates on U(l) 2 \SU(3) arc 



x,y,z € 



0, 



, w, r, t e [0, 2tt] 



(11) 



In this parametrization, the Jarlskog invariant J is given 
by 



1 ■ ■ 2 

J = — sin 2x sin 2z sin y cos y sin w , 

which has a maximal value of — K= w 0.096. 
We obtain a Riemannian measure 

DM := JJ(A - D 3 ) 2 dD l dD 2 dD 3 DU , 



(12) 



(13) 



with 



DU := sin 2x cos 3 y sin y sin 2z dx dy dz dw dr dt , (14) 

on the space of Hermitian 3x3 matrices. We allow these 
to have arbitrary eigenvalues. (In a fermionic mass term, 
the sign of the mass has no physical significance, and only 
m 2 enters in physical quantities.) Note that the subspace 
of matrices with coinciding eigenvalues has measure zero 
and can be ignored in the following discussion. 

Each Hermitian matrix with three distinct eigenval- 
ues is associated with six different elements of M 3 x 
U(1) 2 \SU(3), related by the action of the discrete group 
6 3 : 

U^DU = {U^P- l )PDP~ l {PU) =: U^DU , P G 63 , 

(15) 

where 63 is the symmetric group of degree 3 which per- 
mutes the canonical basis vectors of R 3 . We need to 
consider the space {U(l) 2 x & 3 )\SU(3) x R 3 instead, re- 
stricting the coordinates on U(l) 2 \SU (3) in order to pick 
one of the six matrices related by the 63 action. This 
can be achieved by demanding that the elements of the 
third column have ascending absolute value, 



< y < arctan(sin x) , < x < 



(16) 



Since there are two mass matrices and we need to av- 
erage over two copies of the space of 3 x 3 Hermitian 
matrices, an integral of a quantity such as J 2 becomes 



(J 2 ) = 



J DM DM' e -Tr(M 2 A)-Tr((A/') 2 A') j2( M> M >) 



jDM £>Af'e- Tr ( M2 ' 4 )- Tr (( M ') 2A ') 



(17) 



and J really depends only on V = UU'K We choose J 2 
since all odd powers of J average to zero. We use the 
same parametrization for M' as for M, parametrizing U' 
by x',y',z', etc. 

It should be clear from (fT?]) that multiplying A (or A') 
by a constant is the same as rescaling the eigenvalues 
Di (or D'A, and so amounts to a rescaling of the mass 



3 



scales A and A'. We can therefore, without any loss of 
generality, choose 



where 



.4 



'10 
l/nl 
.0 l//z 2 



A' 



'10 
l/n\ 
.0 



(18) 

where (j, c ,fi u ,fi a and fid are dimensionless parameters 
that we are free to choose so as to reproduce the ob- 
served squared quark masses as expectation values (Df) 
etc. We expect these to be of the same order of mag- 
nitude as the corresponding quark masses, expressed in 
units where A = m% and A' = m&. 

In the case where A and A' are proportional to the 
identity, the exponential factor becomes exp(— aTr(D 2 ) — 
/3Tr((Z?') 2 )), and the integral (fT7|) splits into a product 
of an integral over the eigenvalues, and an integral of J 2 
over ((U(l) 2 x 6 3 )\SU(3)) 2 . The distribution of J 2 is 
then independent of the quark masses, and one would 
expect J 2 to be much larger than observed. A generic 
A couples the diagonal D to U which is an element of 
((U(l) 2 x & 3 )\SU(3)) 2 , and allows the quark masses to 
have an effect on the distribution of J. 



III. RESULTS 

We need to use approximations to evaluate the integral 
P7p. as the expression for J in terms of coordinates on 
{(U(l) 2 x & 3 )\SU(3)) 2 is too complicated to be given 
explicitly. However, since 



it/4 tv/2 



Jd 3 Dd 2 x:= J J{dD k 'Y\{D i -D j f J dx J dz s 2x s 2z 

(22) 

Using the fact that these integrals are dominated by small 
x and z, we find [J] the analytical approximations 



(Dl)*^, {D 2 2 )^^, <^3 2 >«Y' (23) 

and similarly for (D'A 2 . 

We use these approximate analytical results as a con- 
sistency check for more accurate numerical calculations 
of expectation values of squared quark masses. The mass 
scales A and A' are fixed by setting (D 2 ) = (m t /A) 2 and 
((D[) 2 ) — (rrib/A') 2 . By comparing the results obtained 
by numerical integration with the values we want to re- 
produce, we can then fix the parameters fi Cl fi u , fi s and 

Due to the dependence of quark masses on the energy 
scale, described by the renormalization group, there is 
some ambiguity in what is meant by the "quark masses" 
we want to reproduce. Following for example, we 
take all the quark masses evolved to the scale of the Z 
boson mass. We use the central values [6j 



2\ _ K 



(1.27 MeV, 619 MeV, 171.7 GeV) , 
(2.9 MeV, 55 MeV, 2.89 GeV) . (24) 



(m u ,m c ,m t ) 
{m d ,m s ,m b ) 



Comparing our numerical results with (|24|) , we fix the 
parameters to 



Tr(M 2 A) = Tr(D 2 UAU*<) 



D 2 

^1 ■ 2 

-2- sin y 



(19) 



and 



Tv{{M') 2 A') = Tr((D') 2 U'A'(U')^) = ^4 



l\2 

l) -ii 
— sin y 



(20) 

with all remaining terms being non-negative, and ji u <C 
1 and fid <C 1, the integrand is negligibly small unless 
y w and y' w 0. We can therefore do the integrations 
over y and y', using siny « y with the remaining part 
of the integrand taken at y = y' = 0. The result is 
now independent of w and w', and we can integrate the 
expression for J 2 over r,r',t, and t'. Then the integrals 
over both copies of R 3 can be done analytically, but one 
is left with a four-dimensional integral over x,z,x' and 
z' which, in general, can only be evaluated numerically. 

Expectation values for squared masses, however, take 
the relatively simple form 



fd 3 Dd 2 xD 2 ( e -Tr(D*UAUl)\ I 

/r>2\ ~ A / lv=Q (21) 



Jd 3 Dd 2 X (e-T*(D*VAV1) 



Mc = 3 ) Mu = 3 , 

\m t J \m t 

2 3 f m s \ 2 2 12 f m d 

2 V mb J o \m b 



(25) 



to reproduce the corresponding expectation values for 
squared eigenvalues of the mass matrices. (The analytical 
approximation (|23[) is less accurate for the d, s, b quarks 
since the mass hierarchy is milder, so the numerical fac- 
tors given here are numerical fits.) 

We can obtain an analytical expression for the expec- 
tation value of J 2 because the integral ([17]) is dominated 
by the region of small y,y', and z, and therefore to a 
first approximation we only need J 2 at y = y' = z = 0. 
Averaging over r,t,r' and t' gives a factor of 1/2, and 
therefore we can use 



J. 



2 • 2 2-2 / 2 / 2/-4/ 

imall y.y'.z : ~ ~^ Sm x C0S x Sln X COS x C0S z Sm z 



for our calculations. We obtain 

_ 4 

(^small y,y'z) ~ 



Vl5 



m u md mb 



\y=0 



'"■■»*2( 1 + >/f^) 



4 ti t-2 , 



(26) 
(27) 



4 



with 

*i=arctan(^)-^, 
* 2 = a*ctaa(yj:^)-^M, (28) 

where the different prefactors arise from the numerical 
factors given in ([25]) . Note that the top quark mass does 
not appear in this leading approximation. 

Numerical evaluation of (J 2 ) (taken at y = y' = 0, but 
integrated over all z) gives 

AJ = VTJ 2 )" w 7.27 x 10~ 5 . (29) 

Therefore the observed value J « 3.08 x 10 -5 0] is about 
0.42 standard deviations away from zero in our distribu- 
tion, and cannot be viewed as being finely tuned. 

To test the sensitivity of our results to experimental 
uncertainties in the quark masses, we take their values 
at the upper or lower limits given in [|| and obtain the 
corresponding highest and lowest values for (J 2 ). We 
find the bounds 

4.31 x 10~ 5 < AJ < 1.23 x 10" 4 . (30) 

By identifying the dominant regions of the integrals 
appearing in (JT7J) , we are also able to give estimates [J] of 
the expected magnitudes of the components Vtd, Vt s , Vcd 
of the CKM matrix Q , finding 

\V td \ ~Vd~ 0.002, lVi.lw^wO.02, iVdlw — w0.07. 

Ms 

These very rough estimates are within factors of a few of 
the observed values, and reproduce the correct ordering, 
lending further support to the approach taken here. 

IV. CONCLUSIONS AND OUTLOOK 

In this Letter, we have shown that by assuming the 
observed hierarchy in quark masses in a Gaussian distri- 
bution over the space of mass matrices one obtains ex- 
pectation values for J 2 close to the observed value. This 
statistical observation indicates that the mechanism re- 
sponsible the quark mass hierarchy might also explain 
why the observed value for | J| is so small, with no ad- 
ditional fine tuning. In strong contrast, a geometrical 



measure on the space of CKM matrices alone, which does 
not take quark masses into account, makes the observed 
value of J appear highly unlikely [J]. 

These results suggest that the choice of CKM angles 
is closely linked to the observed mass hierarchy. Broadly 
similar conclusions have been reached before [8j , but we 
believe our work represents a significant improvement. 
First, the measure we use is invariant under unitary ro- 
tations of the quark fields Qr and Q' R whereas the mea- 
sure in [8j is not and hence depends on arbitrary con- 
ventions. Second, we have input the observed quark 
masses and obtained fairly tight predictions for (J 2 ), still 
consistent with observation. In contrast, Q assumed 
power law distributions for the magnitudes of the ele- 
ments of M and M' and obtained the much looser pre- 
diction 10~ 8 < |J| < 10~ 4 at one sigma. Finally, we 
have obtained analytical approximations whereas their 
work was entirely numerical. 

Our analysis also applies to the case of massive neutri- 
nos, where the predictions will conceivably be tested by 
future experiments. In the standard theory, the Maki- 
Nakagawa-Sakata matrix [9j which appears is naturally 
an element of U(1) 2 \SU (3), since only phasing of the lep- 
ton charge eigenstates, but not the neutrino mass eigen- 
states (which are assumed Majorana) is possible. Since 
the right phases do not play any role in neutrino oscilla- 
tions, and the relevant J is independent of these phases, 
the calculations are identical to the ones presented here, 
although with the appropriate values of the param- 
eters appearing in A and A'. In the see-saw mecha- 
nism one adds very heavy right-handed neutrinos, and 
the most general mixing matrix would be an element of 
U(1) 5 \SU(6). This is naturally a Kahler manifold, and 
the measure induced by the Kahler metric can be ob- 
tained from the analysis in [lOj. We leave a detailed 
treatment of this case to future work. 

Finally, one could analyze the effects of a fourth gen- 
eration of quarks on CP violation, by repeating the cal- 
culations for 4 x 4 Hcrmitian matrices. If this general- 
ization spoils the agreement with the observed J, one 
might obtain interesting lower bounds on the masses of 
a hypothetical extra generation of quarks. 
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